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The past few years have witnessed tremendous achievements of transformation optics applied to metallic 
plasmonic systems. Due to the poor tunability of metals, however, the ultimate control over surface 
plasmons remains a challenge. Here we propose a new type of graphene plasmonic (GP) metasurfaces by 
shaping the dielectrics underneath monolayer graphene into specific photonic crystals. The radial and axial 
gradient-index (GRIN) lenses are implemented to demonstrate the feasibility and versatility of the proposal. 
It is found that the designed GP-GRIN lenses work perfectly well for focusing, collimating, and guiding the 
GP waves. Especially, they exhibit excellent performances in the THz regime as diverse as ultra-small 
focusing spot (A 0 /60) and broadband electrical tunability. The proposed method offers potential 
opportunities in exploiting active transformational plasmonic elements operating at THz frequencies. 

Transformation optics has been proved to be a general methodology for taming and controlling electromag- 
netic waves in any desired manner 1 . Due to the capabilities to break through the Abbe's diffraction limit and 
implement tight field confinement, surface plasmon polaritons (SPPs) have triggered a large number of 
applications at subwavelength scale 2 " 4 . Transformation plasmonics renders a viable path toward the realization of 
on-chip transformation nanophotonic devices 5 " 8 . Over the past few years, plasmonic gradient-index (GRIN) 
lenses, with the advantages of being flat and free of geometrical aberrations, have been extensively investigated 
theoretically and experimentally 7 " 10 . For instance, Liu et al. theoretically proposed an efficient method to mold 
SPPs bound at metal/dielectric interfaces 7 . They demonstrated that by precisely tailoring the topology of a 
dielectric layer on top of a metal, the propagation of SPPs can be flexibly manipulated. Valle and Longhi reported 
a parabolic GRIN waveguide for subwavelength self-imaging and focusing of the SPP waves, and a cross-connect 
device and an efficient waveguide launcher were numerically demonstrated 9 . Zentgraf et al. experimentally 
validated the feasibility of the plasmonic Luneburg lens and Eaton lens for focusing and bending SPPs by 
modifying the dielectric material adjacent to the metal surface 10 . To date, the majorities of the transformational 
plasmonic devices are on the basis of metal/dielectric model, their limited tunability and increased losses at lower 
frequencies, inevitably, make metallic plasmonics less appealing for applications in the THz and mid-infrared 

1112 

regime ' . 

Since the discovery of easy fabrication processes for remarkably high-quality monolayer graphene by 
Novoselov et al. in 2004, enormous efforts have been rapidly devoted to the emerging field of graphene 13 " 15 . 
With striking electrical and optical properties, a variety of graphene-based devices have been explored, such as 
modulators 16 " 18 , photodetectors 19 " 21 , and cloak 22 . Recently, graphene has been shown to support transverse- 
magnetic (TM) SPP modes at subwavelength scale for frequencies ranging from THz to the visible 23 " 26 . As a 
novel plasmonic material, graphene possesses unprecedented advantages over metals, including flexible tun- 
ability, extremely tight confinement, and low losses at THz and mid-infrared frequencies. Meanwhile, an increas- 
ing number of researches have focused on the hotspot of graphene plasmonics (GPs) 27 " 32 . For example, Vakil and 
Engheta reported that by designing spatially inhomogeneous and nonuniform conductivity patterns, graphene 
could be regarded as a one-atom-thick platform for transformation optical devices 24,25 . Chen et al. experimentally 
investigated the propagation of GP waves in tapered graphene nanostructures using near- field scattering micro- 
scopy 28 . They demonstrated that the wavelength of GP waves was more than 40 times smaller than the wavelength 
of illumination, and the GP waves exhibited extremely small mode volume as well as electrical tunability. All of the 
aforementioned works have theoretically and experimentally proved that graphene can serve as an ideal flatland 
platform for two-dimensional (2D) plasmonic photonic circuits. 
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Figure 1 | Mode index evolution of the guided GP waves, (a) Dependence of GP mode index on the average permittivity of surrounded dielectrics for 
different frequencies with E f = 0.15 eV (corresponds to the gate voltage of 20 V). (b) Dependence of GP mode index on the gate voltage for 
different frequencies with e avg = 4. The Fermi energy above the half of the maximum photon energy (E f > 0.02 eV, corresponds to the gate voltage of 1 V) 
is taken into consideration in (b). The carrier mobility ji = 60000 cm 2 V _1 s _1 is utilized in the calculations. 



In this paper, we propose a flexible design of graphene plasmonic 
quasicrystal (GPQ) metasurfaces by using the specific photonic-crys- 
tal (PC) -like dielectric layer underneath a graphene sheet. Based on 
the proposal, the GP-GRIN lenses acting as the Maxwell's Fisheye 
lens, Luneburg lens, and self-focusing lens are implemented. The 
simulated results show that the designed GP-GRIN lenses work per- 
fectly well for focusing, collimating, and guiding the GP waves. In 
particular, they exhibit excellent performances in the THz regime 
such as broadband electrical tunability and ultra-small focusing spot 
of 1 0 /60 (2 0 is the wavelength of excitation light in vacuum). The 
proposed method allows for new opportunities in exploiting electric- 
ally tunable plasmonic elements for transformation optics at deep 
subwavelength scale. 

Results 

Structure model and analytical theory. We firstly analyze the 
dependence of GP waves on the surrounding environment and the 
Fermi energy of graphene, here considering a single layer graphene 
sheet surrounded by dielectrics above and below graphene with the 
relative permittivities of £ ra and e rb , respectively. According to the 
Maxwell's equations and constitutive relations, the dispersion 
relation for the TM-polarized GP waves propagating along gra- 
phene sheet can be expressed as 23 



trj 0 a g 



(1) 



where p p represents the wave vector of the guided GP waves, p 0 = 
27rMo is me wave vector of excitation light in vacuum, and r] 0 ~377 Q 
denotes the intrinsic impedance of vacuum. The surface conductivity 
of graphene cr g is derived within the random-phase approximation 
and comprised of intraband and interband transitions 33,34 . In our 
analyses, only the intraband contribution is taken into account on 
condition that the half of the photon energy Hco is below the Fermi 
energy of graphene Ef(i.e.,Ef > Hco/2) 35 . Therefore, a g is simplified to 
a Drude form at room temperature (T = 300 K) with Ef ^> k^T 36 



(Ja((D,Ef ,t) = 



ie 2 E { 



(2) 



where co is the angular frequency of excitation light in vacuum, e is 
the electron charge, k B and H are the Boltzmann's constant and 
reduced Planck's constant, respectively. The momentum relaxation 



time t = fiEf/(eVf 2 ) is determined by the carrier mobility fi and Fermi 
energy Ef with the Fermi velocity v f = 10 6 ms" 1 . By substituting Eq. 
(2) into Eq. (1) with fi p >> fi 0 , the propagation constant of the guided 
GP waves is finally obtained as 



3 p (co,E f ,s avg ) = P 0 



2nh 2 (co-\-iT 1 )e a - 
rj 0 e 2 E f 



(3) 



Here e avg = (e ra + s rh )/2 is defined as the average relative permittivity 
of dielectrics above and below the graphene sheet. It can be seen that 
the mode index of the guided GP waves (n p = /? p //? 0 ) is strongly 
dependent on the average relative permittivity e avg and Fermi energy 
Ef besides the excitation frequency co. 

The dependence of GP mode index on the permittivity of sur- 
rounded dielectric is illustrated in Fig. 1(a). It reveals that for a 
certain frequency larger permittivity results in higher mode index 
which is conducive to the achievement of a shorter wavelength of the 
guided GP waves. On the other hand, the dependence of GP mode 
index on the Fermi energy described by the gate voltage is shown in 
Fig. 1(b) (see Methods). It is found that the mode index of GP waves 
can be flexibly tuned from 5 to more than 300 by varying the gate 
voltage. In contrast with metal-based SPPs, the higher mode index of 
GP waves indicates that the deep subwavelength devices and tight 
confinement of energy can be easily realized on graphene. 
Additionally, the propagation loss of the guided GP waves can be 
characterized as Re(/? p ) /Im(/? p ) = cot 23 . In this context, we choose the 
mobility fi = 60000 cm 2 V _1 s _1 and Fermi energy Ef = 0.15 eV 
(corresponds to the gate voltage of 20 V), hence t = 0.9 ps, which 
is rather conservative to feature the practical losses of graphene 37,38 . 
Therefore, for the frequency of 5 THz and the Fermi energy under 
interband threshold, Re(/? p )/Im(/? p ) = 9n, which is quite favorable 
compared with SPPs in metals 23 . 

Based on the aforementioned analyses, we propose a design of the 
GPQ metasurfaces for manipulating the propagation of GP waves on 
the 2D plane. The fundamental structure model of the metasurface 
incorporates a single layer graphene sheet on top of a specific PC-like 
dielectric layer, which consists of periodically arranged two dielec- 
trics with unequal relative permittivities of g rbl (purple) and e r \> 2 
(grey), as presented in Fig. 2(a). A top-gate voltage is employed to 
adjust the carrier concentration in the graphene sheet, which allows 
highly tunability of the GP mode index with great ease by using Eq. 
(5) (see Methods). Owing to the dependence of GP mode index on 
the surroundings, the proposed GPQ metasurface can be equivalent 
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Figure 2 | Fundamental structure model and analyses of the GPQ metasurfaces. (a) Schematic of the GPQ metasurface consisting of four layers: 
substrate, PC-like dielectric, monolayer graphene sheet, and top-gate dielectric, wherein the gate voltage V g is used to adjust the carrier concentration in 
graphene. The PC-like dielectric layer holds two dielectrics with unequal relative permittivities of e r bi (purple) and e r b2 (grey). The excited GP waves 
propagate along the +x-axis direction on the graphene metasurfaces. (b) The analytical (line) and numerical (square) effective refractive indices of the 
GPQ metasurfaces as a function of filling factor at frequencies of 3 (red), 5 (blue), and 7 THz (green) with a gate voltage of 20 V. In the simulation, the 
unit cell size is always kept below one-third of the minimal GP wavelength. The inset of (b) denotes the top view of a unit cell with high (I) and low (II) 
mode index areas, r and a are the radius of the filling area and the period of the unit cell, respectively. 



to a new material with an effective refractive index in accordance 
with the effective medium theory (EMT) 39,40 . The unit cell of the GPQ 
metasurface, as shown in the inset of Fig. 2(b), is divided into two 
sections: high mode index n Y area (I) and relatively low mode index 
riu area (II), and each unit cell can be regarded as a homogeneous 
media within a good approximation. When the excited GP waves are 
guided along the graphene sheet, the effective refractive index of the 
unit cell is governed by 39 ' 40 

«eff =y«i + (i-y)«n. (4) 

Where y = n^/a 2 is the filling factor, r and a are the radius of the 
filling area and the period of the unit cell, respectively. n\ and /% are 
obtained by substituting g avgI = (g ra + e rbl )/2 ande avgII = (e ra + ^hi)l 
2 into Eq. (3). To demonstrate the validity of the EMT approxi- 
mation, the finite -difference time-domain (FDTD) method is 
employed to calculate the effective refractive indices of the metasur- 



faces (see Supplementary Information). According to Eq. (4), the 
effective refractive indices of the homogeneous GPQ metasurfaces 
for different incident frequencies are plotted in Fig. 2(b). It is found 
that the effective refractive indices from FDTD simulations are in 
excellent agreement with the results analytically calculated from 
EMT, which means that the proposed graphene-based structures 
can work as the metasurfaces with the effective indices described 
by Eq. (4). By regularly varying the filling factors of the PC-like 
dielectric layer, consequently, the GPQ metasurfaces with the 
required effective refractive index distributions can be achieved, such 
as the GRIN metasurfaces for transformation optics. 

Numerical results. As proofs of principle, here graded GPQ meta- 
surfaces with the GRIN distributions are implemented by employing 
FDTD method. Two types of classic GRIN lenses are numerically 
investigated on the basis of the graded GPQ metasurfaces: radial and 
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Figure 3 | Radial GP-GRIN metasurface and effective refractive index profiles, (a) Top-view illustration of the radial GP-GRIN metasurface. The filling 
area (purple) gradually decreases along the radial direction, (b) Lateral profiles of effective refractive index of the proposed GP Maxwell's Fisheye and 
Luneburg lenses as a function of filling factor. In the calculations, the frequency of the excitation light and the gate voltage are 5 THz and 20 V, respectively. 
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Figure 4 | 2D effective refractive index maps and numerical results. Left column: the index map (a), field distribution (b), and normalized intensity at 
focusing spot (c) of the GP Fisheye lens. A point source locates on the left rim of the lens. Right column: the index map (d), field distribution (e), and 
normalized intensity at focusing spot (f) of the GP Luneburg lens. A plane wave is launched from the left port of the lens along the +x-axis direction. The 
white dashed circle represents the rim of lens. The FWHMs of focusing spots at x = 5.7 (am along the y-axis direction for the Fisheye lens and 
Luneburg lens are 0.952 and 1.064 (am, respectively. 



axial GP-GRIN lenses. Figure 3(a) schematically depicts the top view 
of the designed radial GP-GRIN metasurfaces. This design can result 
in the achievements of the metasurfaces with the index changing 
gradually along the radial direction. As the examples, 2D 
plasmonic Maxwell's Fisheye and Luneburg lenses based on the 



radial GP-GRIN metasurfaces are demonstrated, which in general 
are difficult to obtain in traditional optical elements 10 . The effective 
refractive index profiles of the graphene-based Fisheye and Luneburg 
lenses, at the frequency of 5 THz, are depicted in Fig. 3(b) according 
to Eqs. (6) and (7) (see Methods). It can be seen that the effective GP 
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Figure 5 | Field distributions of the double-lens system based on two GP 
Luneburg lenses, (a) A plane wave incidents from the left port of the 
system, (b) A point source is set on the rim of the left lens at the position of 
7=0. The black dashed circles denote the rim of lenses. 



mode index varies gradually from 23.6 to 50 as the filling factor 
increases from zero to 0.8. It is worth noting that the index 
distributions of the Fisheye and Luneburg lenses are obtained by 
solely altering the filling factors with other parameters unchanged. 

Being similar to the traditional Fisheye and Luneburg lenses, the 
plasmonic version of the Fisheye (Luneburg) lens can focus a point 
(plane) source to a point image on the rim of the lens 10 ' 40 " 42 . 
Figure 4(a) illustrates the 2D effective refractive index map of the 
GP Fisheye lens. As it can be seen, the background index out of the 
lens equals 23.6 and the diameter of the lens is 11.4 |im. The mode 
index of GP waves at the centre of the lens approaches a maximum 
value of 47.2. In Fig. 4(b), the field distribution shows that when a 
point source is put on the rim of the Fisheye lens the corresponding 
image can be located at the symmetrical position of the rim. For the 
GP Luneburg lens, the background index as well as the lens-size is the 
same as that of the Fisheye lens, besides the maximum mode index of 
GP waves is 33.4 at the centre, as illustrated in Fig. 4(d). Figure 4(e) 
exhibits that a plane wave launched from the left port can be focused 
to a point on the opposite side of the rim. The full-width at half- 
maximums (FWHMs) of the focusing spots for Fisheye lens 
(0.952 |im) and Luneburg lens (1.064 |im), as shown in Figs. 4(c) 
and 4(f), are as small as ~1 0 /60. In fact, at the rim of lens, the effective 
wavelength of GP waves 2 p is ~A 0 /30. The focusing spot size is 
exactly a half of the effective wavelength of GP waves (i.e., A p /2), 
which is the benefit from the high mode index of GP waves. In 
addition, a double-lens system comprised of two Luneburg lenses 
with the rims touched together is simulated. The field distributions in 
Figs. 5(a) and 5(b) prove that the double-lens system can collimate 
the GP waves: only the GP waves propagating along the x-axis dir- 
ection can pass through the system. Furthermore, a four-lens system 
is expected to realize the invisibility of the object in the centre of 
square, which can be regarded as a 2D cloak working similarly as the 
3D cloaks 43 ' 44 (see Supplementary Information). These results indi- 
cate that the GP Fisheye lens and Luneburg lens based on the pro- 
posed approach are potential to function as significant components 



such as collimators, connectors, and cloaks in the deep subwave- 
length photonic circuits. 

Successively, an axial GRIN lens with parabolic refractive index 
distribution is investigated. The axial GP-GRIN lens can be obtained 
by altering the filling factors along the y-axis direction. Figure 6(a) 
illustrates the top view of the axial GP-GRIN metasurface. According 
to Eqs. (4) and (8), the parabolic refractive index profile of the axial 
GP-GRIN lens is plotted in Fig. 6(b) (see Methods). The inset of 
Fig. 6(b) presents the 2D effective refractive index map of the para- 
bolic axial GP-GRIN lens, which is similar to the refractive index 
distribution of the self-focusing fiber 9 . Figure 6(c) shows the self- 
focusing performance of the proposed lens that a plane wave prop- 
agating along the x-direction is gradually focused at the central axis 
(y = 0), and then continuous propagation results in a recovery of the 
plane wave. The focal length of the GP self-focusing lens is j x — 
7.660 |im. Figure 6(d) shows the self-imaging performance that a 
point source located dXy = 2 |im can be transferred to a point image 
dXy = —2 |im, and the focal length is/ 2 — 7.465 |im. The deviation 
between j \ and/ 2 can be attributed to the asymmetric scattering losses 
of the filling areas. The relative deviation ratio, defined as \f x —f 2 \/h p , 
is less than 0.1. Therefore, the deviation is negligible when compared 
with Ap, and the GP self-focusing lens is free of spherical aberration 
for focusing, imaging, and guiding the GP waves. 

Although these lenses are originally designed for the frequency of 
5 THz and the Fermi energy of 0.15 eV, their performances can be 
tuned by employing gate voltage. According to Eqs. (3) and (4), the 
real part of the obtained effective refractive index of the unit cell can 
be denoted as n e f^(D,E { ) = 9co/E { , wherein 0 = 27ih 2 {y(s SLVgL — £ aV gii) 
+ e a vgii}/(^o^ 2 ) is a constant due to that once the lens is constituted 
the lens size and filling factor of each unit cell are kept fixed. On the 
basis of this relationship, n e ^coi,E { )/n e f^(D 2 ,Ef) = coi/co 2 and 
n e ff(<jj,Efi)/n e ff(cD,Ef2) — E^IE^. Thus, the index profile is scaled by 
a factor of a>i/co 2 when the frequency is changed. Similarly, when the 
Fermi energy is tuned, the index profile varies by a factor of E {2 /E {1 . It 
indicates that the index profiles comply with Eqs. (6), (7) and (8) 
throughout the altering of frequency and Fermi energy, which means 
that the proposed lenses are dynamically tunable for wideband fre- 
quencies. Here, we take the GP Luneburg lens for example to dem- 
onstrate this excellent performance. Figure 7(a) depicts the effective 
refractive index profiles of the lens for different gate voltages, which 
agrees well with the above analyses. It is interesting that the GP waves 
are always focused on the rim of lens when the gate voltage is tuned, 
as shown in Figs. 7(b), 7(c), and 7(d), which validates that the 
regained metasurface is still a Luneburg lens. Additionally, it is found 
that when the gate voltage is varied from 15 to 5 V, the FWHMs of 
focusing spots for 4 THz can be reduced from 2.208 to 0.832 |im. 
The conclusion is drawn that the size of focusing spot can be dynam- 
ically tuned by the gate voltage while the centre of focusing spot is 
well maintained on the rim of lens. 

Discussion 

We have proposed and numerically demonstrated a flexible design of 
the GPQ metasurfaces for controlling the SPPs propagating along the 
graphene sheet by shaping dielectric layer underneath monolayer 
graphene into the specific PCs. It is found that the GP-GRIN meta- 
surfaces can be constructed by employing the chirped PC-like dielec- 
tric design. To validate the feasibility and versatility of the proposed 
approach, the graphene metasurfaces with effective refractive index 
profiles of Maxwell's Fisheye lens, Luneburg lens, and self-focusing 
lens are numerically implemented. The simulated results dem- 
onstrate that the designed GP-GRIN lenses exhibit excellent perfor- 
mances in the THz regime such as focusing, collimating, and guiding 
the GP waves. It is found that the spot size is as small as ~X 0 /60, and, 
especially, can be dynamically tuned by adjusting the gate voltage. It 
is interesting that the proposed lenses can work perfectly well for 
broadband frequencies by varying the gate voltage. We also find that 
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Figure 6 | Axial GP-GRIN metasurface and parabolic self-focusing lens, (a) Top-view illustration of the axial GP-GRIN metasurface. The filling area 
gradually decreases along the y-axis direction, (b) Parabolic effective refractive index profile of the proposed self-focusing lens as the function of 
filling factor and width. The inset of (b) shows the 2D effective refractive index map. (c) A plane wave launched from left port is recovered at right port, (d) 
A point source located at y = 2 (am is transferred to a point image at y = —2 (am. 



the lenses can function as significant components such as collimators 
and connectors in the deep sub wavelength photonic circuits. The 
proposed ultrathin planar lenses based on graphene metasurfaces 
are aberration-free and electrically tunable, which are really fascin- 
ating for the active integrated photonic circuits. It is predictable that 
the arbitrary metasurfaces with GRIN index profiles can be projected 
through the proposed method such as optical Janus device 45 , Eaton 
lens 10 , and magnifying lens 46 , which indicates new opportunities in 
developing more active and complex transformational plasmonic 
elements. 

Methods 

Electrical tunability by a top-gate voltage. In the proposed GPQ metasurfaces, the 
carrier concentration in graphene sheet is dynamically controlled by employing 
a top-gate voltage, as shown in Fig. 2(a). The Fermi energy of graphene can be 
estimated using the parallel capacitor model as 47 



\E f (V g )\=Hvf^^V),\ 



(5) 



Where E^V g ) is the Fermi energy at gate voltage of V g and oc c = s 0 s r J{ed) is the 
capacitor constant. In our case, the plate separation of the created capacitor d is 
200 nm and the permittivity of the top-gate dielectric e ra is 3 26 ' 48 . When the gate 
voltage is swept from 1 to 40 V, the Fermi energy experiences a sensitive variation 



covering the region of interest, and hence the GP mode index changes from 5 to more 
than 300, as shown in Fig. 1(b). Therefore, the performance of the proposed GPQ 
metasurfaces can be flexibly tuned by the top-gate voltage through electric field effect. 

Refractive index profiles of the proposed GP-GRIN lenses. The refractive index 
distributions of the traditional Maxwell's Fisheye lens and Luneburg lens are 
expressed as 10,41,42 



n(p) = 



2n 0 



and 



n(p) = n 0 ^2-(p/R) 2 . 



(6) 



(7) 



Where R is the radius of lens, p is the distance to the centre, and n 0 represents the 
refractive index of the background media out of the lens (p > R). The parabolic 
refractive index profile of the self-focusing lens is given by 9 



n(y) = yjn\-g 2 f. 



(8) 



Where g = 2(rii — n 0 2 ) m /w is the gradient constant and w is the full-width size of the 
lens along the v-axis direction. n 0 and n x are the refractive indices at \y\ > w/2 and y = 
0, respectively. For the proposed graphene-based lenses, n 0 is the background 
effective refractive index of the metasurface {i.e., n n ). The filling factor of the jth unit 
cell of the designed lens is calculated by numerically solving the equations n e ff(yj) = 
n(pj) and n^) = n{y^) (j = 0, ±1, ±2, ...). 
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Figure 7 | Performance of the electrical tunability of the GP Luneburg lens, (a) Lateral profiles of effective refractive index of the GP Luneburg 
lens for the gate voltages of 5, 10, and 15 V at frequency of 4 THz. Field distributions and FWHMsofthe focusing spots oftheGP Luneburg lens for 15 (b), 
10 (c), and 5 V (d) at frequency of 4 THz. 



Numerical simulations of the designed GPQ metasurfaces. For the sake of 
simplicity and clarity, the single layer graphene sheet in our model is regarded as an 
ultra-thin metallic film with a thickness of A = 1 nm and an equivalent relative 
permittivity of £ gr) e q = 1 + irj 0 a g /(p 0 A) 24 ' 36 . The relative permittivities of the PC-like 
dielectric layer are are e rb i = 9 and s rb2 = 2 26,49 . In the FDTD simulations, the 
minimum mesh size is set as 0.2 nm, which is small enough to ensure the numerical 
convergence of the results 36 . The unit cell size is kept below one-third of the minimum 
wavelength of the GP waves to prevent the scattering loss. The influences of PC-like 
dielectric on the intrinsic properties of graphene are neglected (see Supplementary 
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